Resonance Tunneling in Double- Well Billiards with a Pointlike Scatterer 
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The coherent tunnehng phenomenon is investigated in rectangular billiards divided into two domains 
by a classically unclimbable potential barrier. We show that by placing a pointlike scatterer inside 
the billiard, we can control the occurrence and the rate of the resonance tunneling. The key role of 
the avoided crossing is stressed. 

KEYWORDS: chaotic tunneling, quantum billiard, delta potential, diabolical degeneracy 



3.65.-W, 4.30.Nk, 5.45.+b, 73.40.Gk 

The tunneling is a quintessentially quantum phe- 
nomenon. However, with the path-integral formulation, 
a systematic semiclassical approximation can be formu- 
lated, and the problem can be treated in terms of the 
instanton, or the "classical orbit with imaginary time" in 
the region classically inaccessible 0. One can assume, 
therefore, that the dynamics of this virtual classical or- 
bits largely determines the characteristics of the tunnel- 
ing phenomena. One can further conjecture that the tun- 
neling in non-integrable system, with its chaotic virtual 
orbit, might display novel features unknown in the text- 
book examples of one-dimensional barrier penetration, 
or more generally, the tunneling in integrable systems 
Recent studies on several model systems show this 
is indeed the case. It has been revealed that very slight 
variation of the shape of the barrier can cause an unex- 
pectedly large enhancement of the tunneling rate [^,^. 
It is also found that the tunneling can be made arbitrary 
small by the proper choice of the shape parameter 
Certain evidence seems to link these characteristics to 
the chaotic character of the virtual classical trajectories 
^^^^ point, however, there is no clear under- 
standing about when and how the large variation in the 
tunneling rate occurs in non-integrable systems, and how 
this might be related to the chaos. 

In this Letter, we start from the premise that we are 
still in need of a model which is rich enough to display 
all the essential features of the non-integrable tunneling, 
yet simple enough to allow the physical intuition to the 
phenomena. As a candidate for such a model, we propose 
a two-dimensional quantum billiard in which a particle 
moves around in a rectangular domain which is divided 
into two sub-domains by a finite-height barrier of rectan- 
gular shape parallel to the outer boundary. We further 
place a pointlike scatterer, or a Dirac's delta potential in 
the domain on and off the barrier. We look at the change 
in the tunneling between two sub-domains while varying 
the location of the delta scatterer. This is intended as a 
"minimal model" for the change of the barrier shape. We 
show in the following that this model enables us to iden- 



tify the key role of the diabolical points and avoided level 
crossing in the non-integrable tunneling. It also reveals 
the condition for the occurrence of the tunneling and its 
enhancement and suppression. The mirror symmetry of 
the system is shown not to be a necessary condition for 
the resonance tunneling. 
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FIG. 1. A schematic depiction of the model billiard sys- 
tem. In the region c — 6/2 < a:: < c-\-h/2, a potential barrier of 
height Uo is placed. The pointlike scatterer located at (X, Y) 
is indicated by a filled circle. 

We consider a quantum particle moving inside a rect- 
angular hard wall of size x Ly whose two sides are 
placed on the x and y axes. We put a barrier potential 
of constant height Uq 



U{r) = Uq9{x2 - x)e[x - xi) 



(1) 



which divides the billiard into the left {x < xi) and the 
right [x > X2) sub-domains. We call the center and the 
width of the barrier c and 6, so that we have xi = c — 6/2 
and X2 = c + 6/2. Wc then place a delta potential of 
strength v a,t r = R. The model is depicted in Fig. 1. 
Formally, the eigenvalue equation is 



-—V^ + U(r)+v5{r-R) 
2m 



tpa{r) = Ealpair). (2) 



It is known that this equation is meaningless as it stands 
pO{ . However, the problem is known to be renornializable 
with the introduction of the formal coupling v in place 
of bare coupling v [^^-|l3| . We write the eigenvalues and 
eigenstates for v = as rjn and which we refer to 
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as unperturbed solutions. The full solution for u ^ is 
given in terms of the unperturbed basis. There are two 
kinds of solutions. The first, a trivial type of solution is 



rjn if 0n(^) = 0, 



(3) 



that is, if R is on the node-line of the wavefunction </<„ (r). 
The second, more generic type of solution of eq. ^ is 
obtained from the equation 



where 



G(f;e) 



G(R;e) - - = 

V 



Vn 



1 



(4) 



(5) 



is the renormalized Green's function. The eigenfunction 
corresponding to the eigenvalue is given by 



E 



MR) 



(6) 



We first consider the case in which the system has mir- 
ror symmetry with respect to the x = line which 
divides the billiard evenly to left and right domains. For 
this, we place the barrier at the center, c = and 
also locate the pointlike scatterer along x — Lx/2 line, 
namely X ~ Lx/2. Because of the mirror symmetry, the 
eigenfunction of the system possesses definite parities. 
When Uq is sufficiently large, there are nearly degenerate 
parity doublets (?Aq-, 4'a+) = (i^a, 4'a+i) with energies 
(Sa-, £a+) = {£a, Ea+i) whcrc wc assume a to be an odd 
integer. From these parity doublets, one can construct 
the wave functions 
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V2 



(7) 



which are localized to opposite sub-domains. The local- 
ized states in each sub-domain Xa^^ evolve according to 



where 



exp {is at) ( x!^^ cos(4^) 



and Aq = £a+ 



(9) 



are the average energy and the energy splitting of the par- 
ity doublet ipa- and ipa+- This means that the system 
oscillates between the states xi^'' and Xa ^ with the half 
period T = tt/A^. This is a quantum beating, or a res- 
onance tunneling between two nearly degenerate states. 
The rate of the tunneling between the left and the right 
domains Da is given by 
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(10) 



Therefore, the energy splitting of the parity doublet Ac 
is a direct measure of the tunneling rate. 
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FIG. 2. The energy splitting of the parity doublets as the 
function of perpendicular location of the pointlike scatterer. 
Its horizontal location is kept at the middle X = Lx/'i. The 
parameters are chosen as = 1.618, Ly = 1/Lx, Uo = 50, 
c = Lx/2, b = Lx/lO and v = 100. 

Because of the two dimensional nature of our model 
system, there is still a freedom to choose the perpendicu- 
lar location of the pointlike scatterer, Y. We look at the 
y-dependence of the energy splitting of the parity dou- 
blets. In Fig. 2, a specific example of Lx = {V5 + l)/2, 
Ly = l/Lx, b = Lx/10 and Uo = 50 is shown. The mass 
of the system is chosen to be m = 27r, and the formal 
coupling is set to v = 100. With this choice of v, the 
states near the ground state come into the strong cou- 
pling region Most notable feature in the Figure is 
the fact that at certain values of Y, the energy splitting, 
thus the tunneling rate is enhanced by several times com- 
pared to the unperturbed case, which, in the Figure, is 
seen as F = and Y — Ly. Another point, no less im- 
gportant, is that there are several values of Y in which 
rhe energy splitting becomes zero. This means that even 
with finite height potential, one can totally suppress the 
tunneling by positioning the pointlike scatterer appro- 
priately. Since the Green's function G{R,e), eq. (5) is a 
monotonously decreasing function of e except at its poles, 
eq. (4) cannot have degenerate solutions. The zero split- 
ting of two energy eigenstates can occur only as the joint 
solution of eqs. (3) and (4). Therefore, the location of 
the pointlike scatterer R* at which the total suppression 
of the resonance tunneling occurs is determined by Q 



(/.„(i?*) = G(i?*,r;„)-i = 0. 

V 



(11) 



Since _R is a two dimensional vector, the solution of eq. 
(11) are the isolated points. These points are known as 
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the diabolical points which play central role in the phe- 
nomena of Berry phase |15| . 
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FIG. 3. First fourteen eigenvalues as the function of hor- 
izontal location of the pointlike scatterer. The perpendicular 
location is fixed at y = 0.4I42Ly. See Fig. 2 for the values 
of parameters. The values X = and correspond to no 
pointlike scatterer. 

Next, we let the horizontal location of the pointlike 
scatterer X to vary. In Fig. 3, we show the first four- 
teen eigenvalues as a function of X while fixing the Y 
at (V2 — One observes that two levels approach 

to each other not only at X = L^/^, but also at other 
places. Indeed, The solution of eq. (11) needs not to 
be restricted to the line x = Lx/2: It can be on any 
nodeline 0„(r) = of the relevant unperturbed wave 
function. Examples are shown in Fig. 4, where we mark 
the diabolical locations with asterisks for several eigen- 
states. When the pointlike scatterer is at the diaboli- 
cal locations, any combination of degenerate states can 
be chosen as the eigenstate. One can form two eigen- 
states which cross over the barrier to have equal amount 
of probability distribution in left and right sub-domains 
with opposite relative signs at one sub-domain. This is 
a generalization of the parity doublet states for the mir- 
ror symmetric system. The doublet states with equally 
partitioned wave function remain to be the eigenfunction 
of the system if the pointlike scatterer is moved off the 
diabolical location in the direction parallel to the y-axis 
along the nodeline 0„(r) — 0. This is understood if one 
notices that one of the unperturbed states remains to be 
the full eigenstate. When the pointlike scatterer is not 
on the nodeline, the eigenstates are localized either to 
the left or the right sub-domains. The numerical exam- 
ple in Fig. 4, which depicts the parameter space (X, y), 
illustrates the situation. When the pointlike scatterer is 
placed at shaded region, the a-th eigenstate is localized 
at the right sub-domain, and when at unshaded region, 
localized at the left. The figure is drawn by defining an 



Iq^ 



dx 



(12) 



with Xd signifying the location of the central nodeline of 
"odd" unperturbed wavcfunction. Typically, Xd is close 
to the center of the barrier, namely, Xd — c. The shade 
indicates the region where la exceeds 1/2. It is on the 
line of I a = 1/2 which separates the shaded and unshaded 
regions where the resonance tunneling can occur. This 
example clearly shows that it is not the mirror symme- 
try of the problem that brings the resonance tunneling 
phenomenon, but rather the existence of the crossover 
between the two wave functions with opposite localiza- 
tion characteristics ll6l . 
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FIG. 4. This chart shows the parameter space {X, Y) with 
the localization status of the a-th eigen function. The shaded 
area represents the location of the pointlike scatterer at which 
the localization index la is more than 1/2. The asterisks are 
the diabolical locations. The resonance tunneling occurs near 
the interface of shaded and unshaded areas. The barrier is 
placed at the middle (c — Lx/2). Other parameters are also 
chosen as in Fig. 2. 

We now consider the more generic case in which the 
mirror symmetry is broken in the unperturbed level. As 
an example, we place the barrier of height Uo = 50 and 
width b = 1/10 X Lx a,t slightly left of the middle c = 
19/40 X Lx- We thus have xi = 17/40 x and X2 = 
21/40 X Lx- Other parameters are kept to be same as 
before. The unperturbed states again form the doublets 
{4'a, 0a-(-i} (ct: odd). But the states and (fia+i are 
localized in the opposite side of the barrier. Therefore, no 
resonance tunneling can be observed in the unperturbed 
system. We proceed to place the pointlike scatterer of 
strength v = 100 at the location {X,Y). In Fig. 5, 
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we show the first fourteen eigenvalues as a function of 
X with a fixed value for Y. {Y = (v^ — l)Ly in this 
example.) One can spot several avoided level crossings. 
Specifically, the ones near X = Lx/'2 appear to separate 
the two distinct regions of the parameter space. 
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FIG. 5. Same as Fig. 3 except tliat the barrier is slightly 
shifted leftward (c = 19/40L^). 
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FIG. 6. Same as Fig. 4 except that the barrier is slightly 
shifted leftward (c = 19/40L^). 

In Fig. 6, the parameter space {X, Y) is depicted with 
diabolical points (asterisks) and localization character- 
istics (shades) in an analogous manner as in Fig. 4. 
Namely, when the pointlike scatterer is placed at shaded 
region, the a-th cigcnstate is localized at tiic right sub- 
domain, and when at unshaded region, localized at the 
left. The line of resonance tunneling on which = 1/2 
now becomes the closed loops which separate the regions 



into left and right localization. The resonance tunneling 
can occur only when R is located very close to this loop. 
When one moves the pointlike scatterer along the loop, 
the rate of the tunneling goes through a variation. The 
rate becomes zero at the diabolical points, and the max- 
imum values tends to come at the points on the loop in 
between the diabolical points in a manner similar to the 
one found in Fig. 2. We conclude that the mirror sym- 
metry is not essential to obtain the resonance tunneling 
in double- well billiards. 

The mechanism behind the localization structure 
shown in Fig. 6 can be further elucidated by a closer in- 
spection of the energy level diagram around the avoided 
crossing region. There, we can assume that considering 
the mixing of two levels is sufficient. For the unperturbed 
doublet {(j)a,<l>a+i} of opposite localization properties, 
we can approximate eqs.(4) and (5) as 
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where we define 

Veff V ^ 



Va+1 



(13) 



2 Va 



2 



- MRY 



(R) 



2 Va+l 



(14) 



Here, we have replaced the energy s in all terms other 
than n = a and n = a -|- 1 by {r]a + 77a+i)/2. The 
quantity v^ff is a function of both a and R. The solution 
of eq. (13) is readily obtained as 



^ _ {Va + Va+l) , Veff{fa + fa+l) 



(15) 



± ^\J{Va+l -r]a+ Vejfifa+1 - /a)}^ + Avlgfafa+l- 

where we have used the notation = (f>a{R)'^. The 
corresponding eigenfunctions are given, apart from the 
normalization, by 

il)oc±{r) oc (j)oc{r) H (t>a+i{r)- (16) 



Along a path of varying R (for example, along the line Y 
= constant as in Fig. 5), the avoided crossing takes place 
at the minimum value of the energy splitting \ea+ — £a- \ ■ 
This occurs when we have 



1 



,{Ry > MiiR) and = 



.(i?)2 



Vejf Va+1 - Va 

for the case of Vgff > 0, and when 

1 _ MijR)^ 

Veff Va+1 - Va 



,{Rf <^ (f>a+i{Rf and 



(17a) 



(17b) 
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for Veff < 0. We consider the case of (a). Analogous 
argument can be made for case (b). Inserting the second 
condition of eq.(17a) into eq. (15), we obtain 



£a± 



where we use the notation 7 = \(j>a+i{R) /4>a{R)\- The 
ratio between the two components of ipa±{r) in eq. (16) 
is given by 



Because we have 7 ~ (the first condition of eq. (17a)), 
this ratio becomes very close to 1, giving one-to-one mix- 
ing of the states of opposite locahzation properties. We 
can conclude that the resonance tunneling occurs at the 
parameter values corresponding to the avoided crossing. 
Therefore, in the whole parameter space {X,Y), the line 
of resonance tunneling is the "valley line" of minima 
of the energy splitting \ea+ — £q-|. At the diabohcal 
points, 7 becomes exactly zero. We also observe from eq. 
(18), that both levels are close to the unperturbed level 
e — T]a+i at the avoided crossing. (This becomes £ = rja 
for Vejf < case.) These features are exactly what one 
finds in the numerical example of Figs. 5 and 6. 

Finally, we discuss our results in a broader view. The 
role of the diabolical points in the resonance tunneling 
has been already recognized by the authors of refs. |^ and 
. The clarification of its distribution in the parameter 
space and its decisive role in shaping the line of resonance 
tunneling is made possible only with the very simple set- 
ting of the current work. Although the model studied 
here is a simplistic one, we can conjecture that our main 
finding is applicable to generic non-integrable tunneling 
systems: The resonance tunneling occurs along the "val- 
ley line" of avoided crossing in the parameter space of 
the barrier shape, and the rate of tunneling undergoes 
a large variation when the shape parameters go through 
the diabolical points. 

In the current study, we have made no immediate ref- 
erence to chaotic motion nor to the semiclassical analysis, 
since the system is solvable in full quantum mechanical 
setting. However, if we look at the doublet states in 
higher energy with higher number of nodes, we expect 
the appearance of complex pattern both in the distribu- 
tion of the diabolical locations and in the shape of the 
line on which resonance tunneling takes place. One can 
also make the tunneling more complex in above sense 
by placing more than one pointlike scatterers in the bil- 
liards. The semiclassical approach to such problem will 
be an interesting and challenging task. 

It is known that the effect of the pointlike scatterer 
disappears in the classical limit |jl^. The fact that the 
classically non-observable object causes the large change 



in the tunneling rate is a good example showing the in- 
tricacy of the quantum physics. We also stress that the 
system considered here is not just of academic nature, 
but may actually be tested either through quantum dot 
( 1 + 7^ ± 7\/4-i-7^) (18) oi' through microwave experiments. Even its application 

as a switching device in the nanometer scale might be 
speculated as a remote possibility. 

In summary, we have studied several characteristics of 
non-integrable tunneling phenomena through the numer- 
ical examples of a double-well billiard with a pointlike 
scatterer. We have identified the key role of the diaboli- 
cal crossings in the tunneling. 

We acknowledge the helpful discussion with the mem- 
bers of the theory division of Institute for Nuclear Study, 
University of Tokyo. 



7± y/A + r 



+ 7 



(19) 
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